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Half-life number sequences collected from nuclear data charts are found 
to obey the Newcomb - Benford law. Based on this fact, it has been sug- 
gested recently, that this law should be used to test the quality of nuclear 
decay models. In this paper we briefly recall how, when and why the 
Newcomb - Benford law can be observed in a set of numbers with a given 
probability distribution. We investigate the special case of nuclear half- 
lives, and show that the law provides no additional clue in understanding 
decay half-lives. Thus, it can play no significant role in testing nuclear 
decay theories. 

PACS numbers: 02.50.Cw, 21.10.-k, 21.10.Tg, 29.85.Fj 

1. Introduction 

The history of the Newcomb - Benford law (NBL) goes back to the nine- 
teenth century, when Newcomb found a rather strange 'law' of nature: the 
end of the logarithmic tables are much less used than their beginning pQ. 
The law had been reinvented later by Frank Benford, who extended it by 
investigating an enormous quantity of data from various sources [2] . To each 
element of a given data sequence he mapped their first significant digit, and 
plotted the histogram of the new sequence. In most of the cases he found 
the same strange pattern: the probability of a number having d as its 
first significant digit follows the law we now call Benford's first digit law: 

P d = \g(l + - d ) (d = l, 2, ...,9). (1) 

We note, that the NBL is more general and covers the other digits, too 
PQ. In spite of this, only the first digit law can be examined in connection 
with nuclear half- lives. The main reason for this is that in many cases we 
do not know the second digit due to experimental errors. What is more, 
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we only know a few thousand half- life values, which does not give enough 
statistics to provide reliable results on the distribution of other digits. 

The NBL was found to be base invariant, meaning that if the numbers in 
the original sequence obeying the law is transformed to a logarithmic base 
of k, and the first significant digits are extracted in this new base, then the 
law still holds: 

P d = \og k (l + -) (d = l, k-1). 

a 

The only obvious limitation is that k cannot be arbitrarily big. 

More excitingly, if one has the elements of a sequence obeying the NBL 
multiplied by a given constant, the new sequence will again follow the law. 
This means that the NBL is scale invariant. 

No matter if the primary data sequence comes from the using of loga- 
rithmic tables or numbers in wealth statistics, magazines, geographic data 
or physics books, there is a very good chance that they will obey the law. 
The validity of the NBL for different kind of number sequences is a long 
standing issue in mathematics and the natural sciences. In physics, there 
are many sequences fulfilling the law approximately well, like physical con- 
stants [3], seismic activity data [1] or the strengths of the lines in atomic 
spectra [5]. The half-lives of a-radioactive nuclei has also been examined 
[6], and most recently this investigation has been extended to j3 decay and 
spontaneous fission, including also ground state and isomeric nuclei 018]. 

In this paper, we examine the recent statements and conclusions regard- 
ing the connection between the NBL and nuclear decay [61 [8] in the light 
of the mathematical explanation of the law. 

2. The chart of half-lives 

2.1. The distribution of the first digits 

In Refs. [3 [71 |8], the authors examine the relative occurrence of the 
first significant decimal digits of nuclear half-life values. The data are found 
to satisfy the NBL. In [6], 477 a decaying nuclei are taken into account, 
while in [8], 2059 j3 decay half-lives are investigated. Ref. [7] examines up 
to 3553 half- lives of nuclei with various decay modes. The situation does 
not change if we add a few hundred or thousand estimated half-life values 
to the number sequence |8] . 

As an illustration, we show the result of our own analysis in Table Q] and 
in Fig. [TJ The input data sequence had been taken from NUBASE2003 |9j. 
The nuclei for which we only know the upper or lower limit of their half-life 
values had been rejected, resulting in an input sequence of 2298 entries. 
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Table 1. The occurrence of the first significant digit of half-lives is in good agree- 
ment with the Newcomb - Bcnford law. 



Digit 


Occurrence 


Expected by the NBL 


1 


701 


692 ± 22 


2 


405 


405 ± 18 


3 


281 


287 ± 16 


4 


210 


223 ± 14 


5 


209 


182 ± 13 


6 


149 


154 ± 12 


7 


112 


133 ± 11 


8 


119 


118 ± 11 


9 


112 


105 ± 10 



} 



} 
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Fig. 1. The distribution of the first significant digit of half-lives follows the New- 
comb - Benford law. 2298 half-life values from NUBASE2003 [S] had been evalu- 
ated. 



The occurrence of the first digits follow the NBL well within two standard 
deviations. 

The calculation of the error of the NBL is based on the binomial distri- 
bution [B] . Let N be the size of the data set (in our case ./V = 2298) and 
Pd the probability of a value having the digit d as its first significant digit. 
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Fig. 2. The result of the ones scaling test executed on our data set. The relative 
number of ones as the first significant digit scatters around its expected value 0.301 
within the error 0.0096. The scaling constant is 1.01. The error is given by Eq.[5J 
and indicated by the horizontal lines. 

One is expected to find N d = NP d values out of N having d as its first 
significant digit, with an error of AfiVj, where 

A[N d ] = ^NP d [l - P d ] (2) 
according to the standard deviation of the binomial distribution. 



2.2. Ones scaling test 

The 'ones scaling test' (OST) is a simple method to test quantitatively 
whether a data set satisfies the Newcomb - Benford law [10j . The data set 
obeys the law if the relative occurrence of the numbers beginning with the 
digit 1 is around Pi = 30.1 %, even after an arbitrary number of multipli- 
cations by a given constant (scaling invariance). The used constant is the 
scaling constant, which is 1.01 in our case. Figure[2] shows the result of 
the OST for the 2298 half-lives we examined. It is clear that the NBL is 
followed quite precisely. 



p08'arXiV printed on June 21, 2010 



■5 



3. The mathematical background of the Newcomb Benford law 

In this section we summarize the conditions when the NBL is satisfied, 
based on the work of Smith [10J. The following reasoning can be used 
whenever one knows the probability density function (p. d. f.) of the set of 
numbers for which the compliance of the NBL is investigated. 

A data set obeys the NBL if and only if its ones scaling test gives lg 2 ~ 
0.301 after any number of multiplications by the scaling constant. Let us 
switch from linear scale to logarithmic one: g(x) = lgx, and let the set of 
numbers have a p. d. f. f(g). The condition that a number has 1 as its first 
significant digit can be described with a sampling function s(g), which is 
a periodic function (a square wave) on the logarithmic scale. Then the Pi 
probability of having 1 as the first significant digit is (Fig.[3|) 

oo 

Pi = J s(g)f(g)dg. (3) 

— oo 

Multiplication by c m (where c > 1 is the scaling constant) works as a 
shift by 7 = m lg c on the logarithmic scale. The area denoted by Pi in Eq.[3] 
remains unchanged if we gain and lose the same area when we perform a 
scaling (Fig. [3]) [TT|. Due to the shape of the sampling function, a scaling by 
10* (i G Z) results in the same Pi, making the results of the OST necessarily 
periodic when c ~ 1 (as in Fig. [2]). According to Newcomb [I], 

Vc G R, c> : 3i G Z, r G [0, 1) : c = 10 i+r . 

Only r plays role in the value of the first significant digit, since i only 
represents the shifting of the decimal point. Newcomb stated that the NBL 
is naturally followed, because the distribution of r is uniform. Though the 
uniformity of r is not at all trivial, by using the train of thought above 
we can see, that his statement can be written in a more precise form: the 
NBL is followed if the distribution of the numbers are log-uniform ranging 
between integer powers of 10. 

To get all functions approximately obeying the law in a mathematically 
rigorous way, the result of the ones scaling test after m scaling iterations 
has to be expressed convolution: 

oo 

Pi (7)= / s(g)f(g- 1 )dg = s(g)*f(-g). 

— oo 

The convolution transforms into a simple multiplication F(/) • S(/) in the 
frequency domain, where F and S are the Fourier transforms of f and s, 
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Fig. 3. Illustration of the sampling of the ones scaling test: a) an example distri- 
bution (log- normal with a mean of 2 and standard deviation 1 on the lg axis), b) 
the result of the multiplication of the number sequence by 10 A9 . 

respectively. Since the support of S(/) is the nonzero integer numbers, and 
the ones scaling test can give a constant 0.301 function only if the F(/) • S(/) 
product is 0.301 at / = and vanishes everywhere else, the following theo- 
rem stands: A set of numbers with the probability distribution i(g) obeys 
the NBL if and only if the Fourier transform F(/) of the probability distri- 
bution function of the numbers vanishes at all nonzero integer frequencies 
(Benford's law compliance theorem, [TO]). For an experimentally good sat- 
isfaction of the law, it is enough for the product to become small at integer 
frequencies. In natural sciences this criterion is most likely to be fulfilled by 
such l{g) distributions the F(/) Fourier transform of which becomes very 
small (F(/) ~ 0) before / = 1, and remains close to zero thereafter. 

4. Conclusions 

The specific distribution of the first digits comes from the method one 
uses when decides whether a number has d as its first digit. This method 
is mathematically defined by the sampling function. It follows, that for any 
given distribution of numbers, the compliance of the NBL (within a given 
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error) can be predicted, but not vice versa: from the compliance of the law 
the distribution cannot be constructed. 

Consequently (and contrary to the findings of [8]), the feeling that 
there is a natural phenomenon behind the NBL is illusoric. While scale 
invariant systems play an important role in today's research, in the case of 
the NBL the source of the scale invariance is the periodicity of the sampling 
function (given the probability distribution fulfills the mentioned criterion) . 

If one knows the distribution of a number set, the NBL carries no addi- 
tional information. Thus, it is only the distribution that should be explained 
or predicted by the theory of a given phenomena obeying the law, and not 
the satisfaction of the law itself. Contrary to the suggestion of e. g. [5j[71[8], 
knowing the p. d. f., the use of the NBL to test a number sequence derived 
from a theory modeling a phenomenon that obeys the law gives no addi- 
tional information on the physics of the system (just as knowing the integral 
of a function on an interval does not help much in finding a function with 
a given shape). 

In the case of nuclear decay, the distribution of the half-lives is known 
(see Fig.|3|). The examined experimental values of the half-lives span 54 dec- 
imal orders of magnitude. It can be seen that this distribution is very close 
to a log-normal distribution, which may suggest a multiplicative process in 
the background. The parameters of the fitted normal distributions can be 
read in Tabled 

By looking at the distribution it is not surprising that half-lives fol- 
low the NBL. The Fourier transform of a Gaussian (see Fig. [5]) with stan- 
dard deviation a is yet another Gaussian with a mean of and a stan- 
dard deviation of l/(27r<r). The standard deviation of the fitted Gaussian is 
3.0139/^2 In 2 ~ 2.56 on the decimal logarithmic scale, making the standard 
deviation of its Fourier transform ~ 0.062. Thus, the transform becomes 
very small at all integer frequencies, resulting only minor fluctuations in the 
value of Pi (7) . 

If a theory of nuclear decay describes this distribution well, it will au- 
tomatically reproduce the Newcomb - Benford law. Thus, the physics is 
behind the form of the distribution, and the models have nothing to do 
with the Newcomb - Benford law itself. 
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Fig. 4. Distribution of the half-lives of the examined 2298 decays. The known half- 
life values cover about 54 orders of magnitude. The binning of the histogram is 
natural logarithm based. 
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Fig. 5. Fourier transform of the half-life distribution and the fitted log-normal 
function. 
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Table 2. The parameters of the Gaussians fitted to the distribution of the half-lives 
of 2298 decays. The parameters differ slightly depending on the resolution of the 
binning. Decimal, natural and binary logarithm based binning has been used. 



binning 


mean 


base mcan 


HWHM 


base HWHM 


lg 


1.0503 ± 0.0051 


11.23 


3.0139 ±0.0060 


1032.52 


In 


3.030 ± 0.018 


20.70 


6.889 ±0.021 


981.42 


lb 


4.609 ± 0.030 


24.40 


9.909 ± 0.036 


961.40 
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